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Abstract — In order to understand the key merits of multiuser diversity 
techniques in relay-assisted cellular multihop networks, this paper 
analyzes the spectral efficiency of opportunistic (i.e., channel-aware) 
scheduling algorithms over a fading multiuser relay channel with K 
users in the asymptotic regime of large (but finite) number of users. 
Using tools from extreme-value theory, we characterize the limiting 
distribution of spectral efficiency focusing on Type I convergence and 
utilize it in investigating the large system behavior of the multiuser relay 
channel as a function of the number of users and physical channel signal- 
to-noise ratios (SNRs). Our analysis results in very accurate formulas 
in the large (but finite) K regime, provides insights on the potential 
performance enhancements from multihop routing and spectrum reuse 
policies in the presence of multiuser diversity gains from opportunistic 
scheduling and helps to identify the regimes and conditions in which 
relay-assisted multiuser communication provides a clear advantage over 
direct multiuser communication. 

I. Introduction 

We consider the uplink and downlink of a cellular multihop/mesh 
system (e.g., IEEE 802. 16j systems), with one base station, one 
fixed relay station and K users. The role of the relay station is 
to enhance end-to-end link quality in terms of capacity, coverage 
and reliability [1], and its presence allows the base station to choose 
between (i) sending/receiving data directly to/from a given user, (ii) 
communicating over a two-hop route where the base station sends 
data to the relay station and the relay station forwards the data to 
the users in downlink, and vice versa for the uplink. We refer to this 
communication model as the multiuser relay channel; which includes 
both the multiaccess relay channel (uplink) [2], [3] and broadcast 
relay channel (downlink) [2], [4]. 

For fixed portable applications, where radio channels are slowly 
varying, multiple access methods based on opportunistic scheduling 
mechanisms take advantage of variations in users' channel quality 
and allocate resources such that the user with the best channel 
quality is served at any given time or frequency (could be subject 
to certain fairness and delay constraints). It has been shown by 
the pioneering works [5] -[7] that the sum capacity under such 
opportunistic scheduling algorithms increases with the number of 
users; yielding multiuser diversity gains by exploiting the time and 
frequency selectivity of wireless channels as well as the independent 
channel variations across users. 

Contributions: While multiuser diversity concepts over traditional 
cellular systems is well understood, there are open issues on the 
design and analysis of resource allocation and opportunistic schedul- 
ing algorithms in relay-assisted cellular multihop networks. Recently, 
low-complexity suboptimal resource allocation algorithms for cellu- 
lar multihop networks were proposed in [8] based on centralized 
scheduling using end-to-end link quality metrics, and were shown 
to simultaneously realize gains from both multiuser diversity and 
multihop relaying to enhance capacity and coverage, provided the 
availability of closed-loop transmission mechanisms. 

In this paper we analyze the spectral efficiency of the opportunistic 
scheduling algorithms (those proposed in [8] as well as others) over 



a multiuser relay channel in the asymptotic regime of large (but 
finite) number of users. In this regime, merely investigating the 
asymptotic capacity scaling as K — > oo (which was the standard 
approach in works on traditional cellular systems) limits the scope of 
the probabilistic analysis and furthermore does not provide accurate 
insights on the large system behavior of the multiuser relay channel 
as the cellular backhaul link (e.g., link between the base station 
and relay station) capacity does not scale with increasing number of 
users. Instead, our approach involves using tools from extreme-value 
theory to approximate the distribution of spectral efficiency, which is 
tight provided that K is large enough. In particular, we characterize 
the limiting distribution of spectral efficiency focusing on Type I 
convergence and utilize it in investigating the large system behavior 
of the multiuser relay channel as a function of the number of users and 
physical channel signal-to-noise ratios (SNRs). Our analysis results 
in very accurate formulas in the large (but finite) K regime, provides 
insights on the role of multihop routing and spectrum reuse policies 
in leveraging the system-level performance of cellular systems in the 
presence of multiuser diversity gains from opportunistic scheduling 
and helps to identify the regimes and conditions in which relay- 
assisted multiuser communication provides a clear advantage over 
direct multiuser communication. 

II. Extreme- Value Theoretic Preliminaries 

Extreme-value theory deals with the stochastic behavior and 
asymptotic distributions of extreme values, such as maxima, min- 
ima or order statistics in general. For instance, the extreme-value 
distributions for maxima are obtained as the limiting distributions of 
greatest values in random samples of increasing size, and are used as 
an approximation to model the maxima of long (finite) sequences of 
random variables. To obtain a non-degenerate limiting distribution, it 
is necessary to reduce the actual greatest value by applying a linear 
transformation (or any transformation) with normalizing coefficients 
which depend on the sample size. 

Let £1,^2, £m be independently and identically distributed 
(i.i.d.) random variables drawn from a common cumulative distribu- 
tion function (CDF) F(x) and denote the maximum of the sequence 
by Xm ~ niax m= i i ... i ]\f £ TO . If there exist sequences of constants 
am > 0, bjwr, and some nondegenerate distribution function /i such 
that (Xm — &m ) /a m converges in distribution to fi as M — * oo , then 
fi belongs to one of the three families of extreme-value distributions 
[9]: Frechet, Weibull and Gumbel distributions. The distribution 
function of £ m , F, determines the exact limiting distribution. Hence, 
in the asymptotic regime of large M, it is not necessary to know the 
detailed nature of F in order to characterize the statistical behavior of 
Xm', and it is sufficient to know which limiting distribution (if any) 
F gives rise to, i.e. to which one of the three domains of attraction 
F belongs. 

In [10]-[11], von Mises (1936) and Gnedenko (1943) determine 
the necessary and sufficient conditions for weak convergence of Xm 



to each of the three types of asymptotic limiting distributions, given 
the parent CDF F corresponding to £ m . These conditions involve the 
tail behavior of 1 — F(x) as x increases - for each type of limit. For 
the case in which F(x) has a density function f(x); the sufficient 
conditions are due to von Mises [10]. 

In this paper, we shall mainly be concerned with the case for 
which wireless fading results in received channel power distributions 
following the Type I extreme-value distribution, where fj.(x) = 
exp(— exp(— x)); and this distribution function is known as the 
Gumbel distribution. Assuming absolutely continuous parent distri- 
butions F with density /, and that there exists a real number x\ 
and X2 < oo with X2 = sup{x : F(x) < 1} such that, for all 
xi < x < X2, f has a negative derivative /' and f(x) = for 
x > X2, we can use the following sufficient condition to determine if 
the parent distribution function F(x) belongs to the Type I domain 
of attraction (Theorem 1.6.1 in [9]): 
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Many fading distributions, e.g. Rayleigh, Rician, lognormal, could be 
given as examples leading to Type I convergence. We can determine 
the normalizing constants ajw and 6m for Type I convergence, by 
solving for &m in 1 — F(Pm) = l/M and setting cim = ??(6m), 
where r\ is the reciprocal hazard function r\{x) = (1 — F(x))/ f(x). 

Defining the spectral efficiency function as C(x) — log 2 (l + x), 
the following lemma extends the convergence results on maxima 
of {£m}m=i to tne maxima of spectral efficiencies by determining 
the constants necessary to preserve the same limiting extreme-value 
distribution. The simple proof (omitted) exploits the fact that C(x) 
is a monotonically increasing function of x. 

Lemma 1: Assume that there exist sequences of constants a,M, 6m 
(the choice will depend on the distribution F) such that 
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as M —* oo, for some limiting extreme-value distribution fj,. Now, 
let Ym = max m C(SNR£ m ) for a given constant SNR > 0. If (0 
holds, then there exist sequences of constants Cm, dm such that 



dt. 



CM 



where 



CM 



log 2 (e)SNRa M 



I + SNR&m 
provided that cm — > as M — * oo 



< x 



p,{x) 



log 2 (l + SNR bi, 



(3) 



III. Multiuser Relay Channel Model 

Consider the network depicted in Fig. [T] with K + 2 nodes, in 
which K users, indexed by k — 1,...,K, send/receive information 
to/from a base station. The relay station is designated to help users 
transmit/receive information utilizing its high capacity backhaul link 
to the base station. In other words, multiple users share a single relay 
for uplink and downlink (i.e. multiaccess and broadcast). 

We assume that all the links over the multiuser relay channel are 
corrupted by additive white Gaussian noise (AWGN). Furthermore, 
the links between the base station and users are assumed to be under 
frequency-flat multiplicative fading i.i.d. across users, with complex 
channel gains {hk}£ = i, where hk G C is a complex random variable 
representing the channel between the base station and user k, drawn 
from an arbitrary continuous distribution Fh with E[|/ifc| 2 ] — 1, Vfc. 
The average received SNR between the base station and user k 
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Fig. 1. Multiuser relay channel model (both uplink and downlink). 



equals SNR^'. Analogously, the links between the relay station and 
users have average SNRs given by {SNR^, r ^}^ =1 [3 and are under 
frequency-flat fading i.i.d. across users, with complex channel gains 
{gk}k=i> where gu £ C represents the channel between the relay 
station and user k, drawn from an arbitrary continuous distribution 
F g with E[| cjffe | 2 ] = 1, Vfc. The set of channels {hk}h=i an d {SfcjfcLi 
are independent. It is assumed that the cellular backhaul link between 
the base station and relay station is an AWGN line-of-sight (LOS) 
connection with received SNR equal to SNRs. The fading states 
over the multiuser links remain constant during the transmission 
of a codeword (slow fading assumption). Due to slow fading, each 
terminal in the multiuser relay channel is considered to obtain full 
channel state information regarding its transmission/reception links. 

For the following analysis, we assume a time-division based (half 
duplex) relaying constraint for multi-hop routing protocols, which is 
due to the practical limitation that terminals can often not transmit 
and receive at the same time. In particular, we consider a two-phased 
decode-and-forward protocol; where, for a given routing path between 
the base station and a given user, the relay station hears and fully 
decodes the transmitted data signal in the first phase and forwards 
its re-encoded version in the second phase. 

IV. Relaying in the Presence of Multiuser Diversity 

Consider the scheduling problem such that K users in the multiuser 
relay channel are to be assigned time-slots for transmission/reception 
over a common bandwidth. In particular, we are interested in the 
max-route centralized scheduling algorithm proposed by [8], which is 
an extension of the maximum signal-to-interference-plus-noise ratio 
(max-SINR) algorithm to relay-assisted networks. Analogous to the 
max-SINR scheduler, which always serves the user with the best 
instantaneous rate at any given time/frequency resource, the max- 
route scheduling algorithm schedules the user with the highest end- 
to-end instantaneous rate. More specifically, the algorithm works as 
follows in the multiuser relay channel: (i) The base station compares 
the instantaneous rates over the direct links between itself and the 
users determined by {ftfe} fc=1 and finds the best user, (ii) The relay 
station compares the instantaneous rates over the links between itself 
and the users determined by {gk}k=i an d informs the base station 
about the index and instantaneous rate corresponding to the best user, 
(iii) The base station compares the link quality of the best user over 
its direct link with the end-to-end link quality of the best user over 
the relay-assisted multihop route (accounting for the transmission 
over the cellular backhaul connection between the base station and 
relay station assuming equal time-sharing between relay-user link 

'Superscript b stands for "base" station and r stands for "relay" station. 



and wireless backhaul link), and schedules the user with the highest 
instantaneous rate. Here, we have dropped the dependence of average 
received SNRs on the user index k and set SNR^' = SNR (6) and 
SNR<: r) = SNRW, Vfc. 

Assuming Gaussian inputs, i.e., all input signals have the tempo- 
rally i.i.d. zero-mean circularly symmetric complex Gaussian distri- 
bution, the maximum supportable end-to-end spectral efficiency over 
the multiuser relay channel achieved by the max-route opportunistic 
scheduling algorithm is given by (in bits per second per Hertz 
(bps/Hz)) 



= max |maxC(SNR w |/i fc r), 
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where the 1/2 penalty factor for the end-to-end spectral efficiency 
of the relay-assisted multihop route from the base station to the user 
is due to the half-duplex constraint imposed at the relay station. In 
10, we have two maxima among K i.i.d. spectral efficiency random 
variables, however the maximum corresponding to the relay-user link 
is truncated due to the presence of the constant capacity backhaul link 
between the base station and relay station. The end-to-end spectral 
efficiency is represented as the maximum of the spectral efficiency 
of the direct link and that of the relay-assisted two-hop link, which 
implies that the users are scheduled in the presence of a channel- 
adaptive (i.e., dynamic) multihop routing mechanism which finds 
their best route to the base station (i.e. switching between direct 
transmission and two-hop relaying). 

As the result of Lemma 1, in the asymptotic regime of large K, 
there exist sequences of constants c^, n ^ 
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where the independent random variables Q^ h ' and 0' 9 ' follow the 
extreme-value distributions pS h ' and yp e ' (which may be the same 
or different), respectively, which in turn are determined by fading 
distributions Fh and F g . 

Theorem 1: In the asymptotic regime of large K, assuming that 
Fh and F B both belong to Type I domain of attraction (i.e. both 



fjP 1 ' and /i y91 follow the Gumbel distribution), the probability that 
the direct link yields higher spectral efficiency than relay-assisted 
two-hop routing is specified by Pc = liniK^oo Vk such that 
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where events A and C are defined as 
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the corresponding probabilities are given by V(A) = e ZK < 2 t 
P (C|.4 C ) = 1 — e _Zjf ' 3 and equation ©, with sequences zk,i, zk.2 

2 -i-> denotes convergence in distribution. 
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Fig. 2. Probability that the direct link yields higher spectral efficiency than 
relay-assisted two-hop routing as a function of the number of users. 



and zk 3 such that 
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and function Q defined as Q(x) — e v ^ dy. 

Assuming Rayleigh fading distribution on Fh and F g , and setting 
SNR (b) = 1, SNR (r) = 10 and SNR S = 100, we plot in Fig. [2] 
empirically generated (based on Monte Carlo simulations) probability 
values for as a function of the number of users K. We also 
plot our analytical result in Theorem 1 obtained by using extreme- 
value theory. We observe that our analytical results are well in 
agreement with the empirical results, and that higher level of accuracy 
is achieved with higher K. Moreover, this comparison reveals an 
interesting tradeoff on the merits of multihop relaying; showing that 
relaying becomes less beneficial as more multiuser diversity gains can 
be realized. In other words, the direct link spectral efficiency grows 
unboundedly with increasing number of users while the end-to-end 
spectral efficiency of the relay-assisted two-hop route converges to 
C(SNR_g)/2 due to the presence of the constant capacity backhaul 
link, designating the use of the direct link more advantageous as K 
grows larger. However, in many practical deployment scenarios we 
have SNRs > SNR (,,) and SNR S > SNR (r) ensuring that the 
wireless backhaul links have high enough capacity so that they do 
not become the bottleneck for the multihop communication in the 
presence of opportunistic scheduling; enabling simultaneous realiza- 
tion of multiuser diversity gains along with capacity and coverage 
enhancements from relaying techniques particularly for cell edge 
users for which SNR (r) > SNR (b) . 

V. Spectrum Reuse in Multiuser Relay Channels 

In this section, with cellular applications in mind, we will focus 
on a special mode of relay-assisted downlink communication over 
the multiuser relay channel depicted in Fig. [3] with U + V + 2 nodes 
(K = U + V users). The users are divided into two categories: U 
far users, indexed by u — 1, ...,[/, with poor quality links to the 
base station (e.g., users at cell edge) and V near users, indexed by 
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Fig. 3. Multihop broadcast channel model (downlink). 
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Fig. 4. Multihop time allocation policies based on orthogonal and simulta- 
neous transmissions. 



v = 1, V, with high quality links to the base statiotjj The role of 
the relay station is to enhance end-to-end link quality for the far users 
in terms of capacity and coverage using multihop routing techniques 
[1] and its presence allows the base station to communicate with the 
far users over a two-hop route where the base station sends data to the 
relay station over a high capacity wireless backhaul link and the relay 
station forwards the data to the far users with possible interference 
from the base station over the second hop. Meanwhile, the near 
users receive downlink data from the base station directly (with no 
help from the relay station) over the same bandwidth with possible 
interference from the relay station. We refer to this communication 
model as the multihop broadcast channel; which is a special case of 
the more general broadcast relay channel. 

While we adopt all the general assumptions in Section II on 
channel modeling and statistics, we now introduce a more specific 
setup to simplify analysis and address the distinction between far 
and near user models. In particular, for the links between the base 
station and users, we drop the dependence of average received SNRs 
on the user index k by setting SNR^ 6 ' = SNRj? for U far users 
and SNRl b) 



k ~ -"""F 
SNR^ J for V near users, and similarly for the 



(r) 



snrP 



links between the relay station and users, we set SNR 
for U far users and SNR^ r ' = SNR^' for V near userfl For 
notational clarity, we re- write the complex channel gains as follows: 
Channel links between the base station and far users are represented 
by {/if,u}u=i, those between the base station and near users are 
represented by {hN,v}v-i, those between the relay station and far 
users are represented by {jf,u}u=i, and those between the relay 
station and near users are represented by {<?jv,„}Jf = i. 

Consider the scheduling problem such that U far users and V near 
users are to be assigned time slots for reception over a common 
bandwidth. This problem involves transmissions over three types of 
links: (i) Lb'- Wireless backhaul link between the base station and 
relay station, (ii) L f '■ The link between the relay station and far users 
and (iii) Lm'- The link between the base station and near users. We 
assign positive time-sharing coefficients /3b, Pf and Pn to links Lb, 
Lf and Lm, respectively, such that Pb + Pf + Pn = 1 and define 
the following time-allocation policies as depicted in Fig. [4] 

3 In contrast to Section IV, this section assumes that a network entry and 
handoff algorithm has been executed prior to the designation of far and near 
users and that there is no dynamic routing mechanism for a given user to 
switch between direct link and relay-assisted two-hop route. 

4 Subscript F stands for "far" and N stands for "near". 



a ) Orthogonal transmission (no spectrum reuse): Links Lb, Lf 
and Ljv are active over different time resources with the corre- 
sponding time-sharing constants [3b, Pf and Pn- The relay station 
compares the channel gains {<?f,u}u=i °f me far users and assigns 
link Lf for downlink transmission to the far user with the highest in- 
stantaneous rate. Analogously, the base station compares the channel 
gains {h,N,v}X=i of the near users and assigns link Ljv for downlink 
transmission to the near user with the highest instantaneous rate. 

h) Simultaneous transmission (spectrum reuse): Link Lb is ac- 
tive over Pb fraction of the time, while links Lf and Ln are 
simultaneously active over Pf + Pn fraction of the time. The relay 
station accounts for average received signal-to-noise ratios SNR^\ 
SNR^ and channel gains {</f,u}u=i, {h-F,u}u=i to schedule link 
Lf for downlink transmission to the far user with the highest 
instantaneous rate. Similarly, the base station accounts for SNR^', 
SNR^', {h,N,v}X=ii {gN,v}X=i t° schedule link Ljv for downlink 
transmission to the near user with the highest instantaneous rate. No 
coordination is assumed to be present between the base station and 
relay station to manage the resulting intracell interference. 

We now study the spectral efficiency performance of these two 
downlink transmission protocols over the multihop broadcast channel 
in the asymptotic regime of large (but finite) U and V: 

A. Orthogonal Transmissions 

The maximum supportable end-to-end spectral efficiency over the 
multihop broadcast channel achieved by the max-SINR opportunistic 
scheduling algorithm in the presence of the orthogonal transmission 
protocol is given by (in bits per second per Hertz (bps/Hz)) 
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Our goal is to characterize the extreme-value distribution of C ort 
in the limit of large U and V, i.e., large number of far/near users. 
As the result of Lemma 1, in the asymptotic regime of large U, V, 
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Theorem 2: In the asymptotic regime of large U, V, assuming 
that Fh and F g both belong to Type I domain of attraction (i.e. both 
pP 1 ' and // 9 ' follow the Gumbel distribution), the average spectral 
efficiency of the orthogonal transmission protocol defined as H ort = 
E [C ort ] can be expressed as 



H° rt =[3 N (c^K + d^+PB C(SHR B )-f3 F 4?' EiO,: 
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where k ~ 0.57721566 is Eider's constant (since E[0 
under Type I convergence), Ei(x) is the exponential integral function 
defined by Ei(rz;) = ^— dy and 
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B. Simultaneous Transmissions 

The maximum supportable end-to-end spectral efficiency over the 
multihop broadcast channel achieved by the max-SINR opportunistic 
scheduling algorithm in the presence of the simultaneous transmission 
protocol is given by (in bits per second per Hertz (bps/Hz)) 
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where SINRf, u is the SINR for far user u and SINRjv^ is the SINR 
for near user v given by 
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As the result of Lemma 1, in the asymptotic regime of large U, V, 
there exist sequences of constants p\j\ o\j\ Py, such that the 
limiting distribution of C slm is expressed as 

C sim (13 F + p N ) ( P i n) e (n) + 4 n) ) 

+ min \p B C(SNR s ), (0 f + p N ) (^6« + 

where random variables G^' and 0'"' follow the extreme-value dis- 
tributions and fv- n \ respectively, resulting from the convergence 
of the maxima of {SINRf,u}u=i and {SIN Rjv,d }Jf=i - 

Theorem 3: In the asymptotic regime of large U, V, assuming 
that both (J,") and pS n ' follow the Gumbel distribution, the average 
spectral efficiency of the simultaneous transmission protocol defined 
as S slm = E [C slm ] can be expressed as 
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Assuming Rayleigh fading distribution on Fh and F g and setting 
SNR S = 1000, SNR^' = SNR^ = 100, SNR^ 6) = SNR^ = 1, 
we plot in Fig.|5]average spectral efficiency as a function of the num- 
ber of far / near users (set U = V, K = 2U = 2V) for orthogonal 
and simultaneous transmission protocols over the multihop broadcast 
channel with time-sharing coefficients set to p B — 0.25, Pf = 




Number of far/near users 



Fig. 5. Average spectral efficiency as a function of the far / near users for 
the orthogonal and simultaneous transmission protocols. 



0.25, Pn = 0.5. Here, we compare empirically generated average 
spectral efficiencies H olt and H slm with their analytical counterparts 
in J 10b and d 1 21 . The empirical results are obtained by averaging 
the expressions in l|9} and Jilt over a large number of randomly 
generated fading realizations (based on Monte Carlo simulations). 
From Fig. [5] we validate that our analytical results are well in 
agreement with the empirical results and that higher level of accuracy 
is achieved with higher U, V. Moreover, we observe that both orthog- 
onal transmission and simultaneous transmission protocols realize 
multiuser diversity gains from opportunistic scheduling techniques 
due to the increase of the average spectral efficiency in the number 
of far / near users. Finally, we find that significant spectral efficiency 
gains can be achieved through spectrum reuse by the base station and 
relay station based on the superior performance of the simultaneous 
transmission protocol over the orthogonal transmission protocol. We 
emphasize that such spectral efficiency gains are possible despite 
that the simultaneous transmission protocol requires no scheduling 
coordination between the base station and relay station and does not 
enforce any policy for interference management. 

References 

[1] O. Oyman, J. N. Laneman, and S. Sandhu, "Multihop relaying for 
broadband wireless mesh networks: From theory to practice," IEEE 
Communications Magazine, 2007, to appear. 

G. Kramer, M. Gastpar, and P. Gupta, "Cooperative strategies and 
capacity theorems for relay networks," IEEE Trans. Inf. Theory, vol. 



[2] 



51, no. 9, pp. 3037-3063, Sep. 2005. 

nd J. N. Laneman, ' The divs 



[3] D. Chen arid J. N. Laneman, "The diversity-multiplexing tradeoff for the 
multiaccess relay channel," in Proc. Conference on Information Sciences 
and Systems (CISS), Princeton, NJ, March 2006. 
[4] Y. Liang and V. V. Veeravalli, "Cooperative relay broadcast channels," 

IEEE Trans. Inf. Theory, 2007, to appear. 
[5] R. Knopp and P. Humblet, "Information capacity and power control in 
single cell multiuser communications," in Proc. IEEE Int. Computer 
Corf,, Seattle, WA, Jun. 1995. 
[6] D. Tse and S. Hanly, "Multi-access fading channels: Polymatroid 
structure, optimal resource allocation and throughput capacities," IEEE 
Trans. Inf. Theory, vol. 44, no. 7, pp. 2796-2815, Nov. 1998. 
[7] P. Viswanath, D. N. C. Tse, and R. Laroia, "Opportunistic beamforming 
using dumb antennas," IEEE Trans. Inf. Theory, vol. 48, no. 6, pp. 
1277-1294, Jun. 2002. 
[8] O. Oyman, "OFDM 2 A: A centralized resource allocation policy for 
cellular multi-hop networks," in Proc. IEEE Asilomar Conference on 
Signals, Systems and Computers, Monterey, CA, Oct. 2006. 
[9] M. R. Leadbetter, G. Lindgren, and H. Rootzen, Extremes and Related 
Properties of Random Sequences and Processes, Springer- Verlag, New 
York, NY, 1983. 

[10] R. von Mises, "La distribution de la plus grande de n valuers," 
Revue Math, de 1'Union Interhalkanique, vol. 1, pp. 141-160, 1936, 
Reproduced in "Selected papers of Richard von Mises", Amer. Math. 
Soc, vol. 2, pp. 271-294, 1954. 

[11] B. V. Gnedenko, "On the limiting distribution of the maximum term 
in a random series," Breakthroughs in Statistics, vol. 1, pp. 195-225, 
1992, Springer- Verlag, edited by S. Kotz and N. L. Johnson (original 
paper published in 1943). 



